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ABSTRACT 

We present a comparison of archival Spitzer spectra of the strongly variable T Tauri EX Lupi, 
observed before and during its 2008 outburst. We analyze the mid- infrared emission from gas-phase 
molecules thought to originate in a circumstellar disk. In quiescence the emission shows a forest of H2O 
lines, highly excited OH lines, and the Q branches of the organics C2H2, HCN, and CO2, similar to the 
emission observed toward several T Tauri systems. The outburst emission shows instead remarkable 
changes: H2O and OH line fluxes increase, new OH, H2, and Hi transitions are detected, and organics 
are no longer seen. We adopt a simple model of a single-temperature slab of gas in local thermal 
equilibrium, a common approach for molecular analyses of Spitzer spectra, and derive the excitation 
temperature, column density, and emitting area of H2O and OH. We show how model results strongly 
depend on the selection of emission lines fitted, and that spectrally-resolved observations are essential 
for a correct interpretation of the molecular emission from disks, particularly in the case of water. 
Using H2O lines that can be approximated as thermalized to a single temperature, our results are 
consistent with a column density decrease in outburst while the emitting area of warm gas increases. 
A rotation diagram analysis suggests that the OH emission can be explained with two temperature 
components, which remarkably increase in column density in outburst. The relative change of H2O 
and OH emission suggests a key role for UV radiation in the disk surface chemistry. 

Subject headings: circumstellar matter — molecular processes — stars: activity — stars: individual: 
(EX Lupi) — stars: pre-main sequence — stars: variables: T Tauri 



1. INTRODUCTION 

EX Lupi is a young star-|-disk system whose photo- 
me tric variability was discovered in 1944 (|McLaughlin| 
1946). Since then the source has been monitored and 
showed repetitive eruptive phenomena related to signif- 
icant ch anges in the amount of accreting material onto 
the star ( Lehmann et al.|1995 Herbig|2007 Juhasz et al. 



20 1 1 1 . Spectroscopic observations during the quiescent 



phases between the recurring outbursts suggest s imilar- 



ity with classical T Ta uri stars of MO type (Herbig et al. 
2001 Sipos et al.|2009 ). The discovery of sirnilar variable 



is just the most recent of the series of eruptive events 
recorded during the half- century- long monitoring of the 
source, but more interestingly it is the most extreme ever. 
It showed an increase of ^5 mag in visual light over seven 
months (January-September 2008), compared to the typ- 
ical increase of ~l- 2 mag observed in its previous char- 
acteristic outburst s (|Aspin et al. [ 2010 1 as well as in other 
T Tauri stars (e.g. Herbst eraTjlUMl r A comparable op- 
tical brightening of the source was previously observed 
only in 19 55, but spect roscopic observations were not 



sources conhrmed EX Lupi as the prototype of a class 
of young strongly a ctive T Tauri stars, called "EXors" 



(Herbig 1989 20081. This class of objects has been pro- 



posed as an mterm ediate stage between FUors and clas 



sical T Tauri stars (Herbig"1977 



exhibiting repetitive outbursts t 



Teodorani et al ||l999[ ) 



of shorter duration than FUors ( |Herbig|2007| [2008 ) . The 



2008 outburst of EX Lupi, that we stiidy in this work 
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performed ('Herbigl|1977|. The 2008 event, instead, was 
monitored at several wavelengths over its entire duration, 
which has provided a deeper insight into the complex 
behaviour of EX Lupi and the still not wel l understood 
me chani s m producing its strong outbursts ( Abraham et 
al. 2009: 'Aspin et al.''2010'; 'Gro sso et aL]|2(jiO{ [Goto^ 

In the last tew years tEe Spitzer Space Telescope 
( Werner et al.|2004 ) has enabled the study of the molec- 
ular content in circumstellar disks through their mid- 
infrared (MIR) emission, unveiling a rich chemistry in 



lat are both weaker and al.,2011, Juhasz et al.,2011^; Kospal et al.||2011 1 . 



2008 



the innermost re gions of T Tauri d isks ( Salyk et al 

f^OlT; 'Carr & Najita 2008, 20lT| |Pascucci ei al.||2009 
|Pontoppidan et al.^^2010a) . So far, such emission has 
been fo und to be common in the T Tauri systems ob- 
served ( Pontoppidan et ar]|2010a ) and has been studied 
by means of simple single-slab models that assume the 
emitting gas to be in local thermal equilibrium (LTE). A 
certain diversity in the molecular emission ob served to- 
ward T Tauri systems is apparent ([Pa scucci e t al. | |2009 



Pontoppidan et al. 2010a Carr & Najita ,2011 ), anc 



disks in more evolved systems as well as around Her 
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TABLE 1 

Summary of EX Lupi Data from the Spitzer Archive (IRS Staring Mode). 



Obs. ID 


Date 


Phase 


Investigator (identifier) 


Module 


Ramp 


tint (s) 


1 79 




c(iiieseciicc 




on 
LH 

SL2, SLl, LLl 


•^n V 1 

OU X 1 

60 X 1 
14 X 1 


fin 
120 
24 


3716 


2005-03-18 


quiescence 


Stringfellow (S05) 


SH 
LH 

SL2, SLl 


120 X 2 
60 X 2 
14 X 4 


480 
240 
112 


477 


2008-04-21 


outburst 


Abraham (AOS) 


SH, LH, SL2, SLl 
SH, LH (backgr.) 


6x4 
6x4 


48 
48 


50641 


2008-05-02 


outburst 


Carr (C08) 


SH 
LH 

SH (backgr.) 
LH (backgr.) 


30 X 4 
60 X 4 
30 X 2 
60 X 2 


240 
480 
120 
240 


524 


2008- 10-10 

2009- 04-07 


quiescence 


Abraham (A09) 


SL2, SLl, LL2, LLl 
SL2, SLl, LL2, LLl 


6x4 
6x4 


48 
48 



big AeBe stars seem to be depleted in warm molecules 



compared to T Tauri disks (Pontoppidan et al. 2010a 
Najita et a"Ll|2010| [Fedele et al.||2011| ). J^'rom the theo 



retical side, several attempts have been made to explain 
the formation and survival of molecules (especially wa- 
ter) in disks and their abundance during the early phases 



et al.|2009| 


Bethell & Bergin,2009j 


Vasyunin et al.|2011pT' 


How does t 


le warm molecular gas composition depend on 



stellar -f-disk properties, and how may the interplay and 
evolution of the different components enable (instead of 
hinder) the availability of essential molecules for planet 
formation? The most recent outburst in EX Lupi has 
provided an extraordinary opportunity to address this 
compelling questions. While it is difficult to identify one 
factor that is the main responsible for the diversity in 
molecular emission of classical T Tauri systems, due to 
the many properties that vary from system to system, 
in the case of EX Lupi we can study a single system 
where only one parameter (the stellar -|-accretion lumi- 
nosity) dramatically increases, and in a timescale that 
is feasible for us to monitor. Crystal formation in out- 
bu rst from amorphous du st particles was recently shown 

by ■ 



Abraham et al. 



( 2009 1 comparing Spitzer spectra of 

EX Lupi, suggesting that episodic increases in disk tem- 
perature during outbursts may contribute significantly to 
the formation of some ingredients that are now present 
in our Solar System. Using the same spectra, but fo- 
cussing on the gas, we show here how H2O (water) emis- 
sion varies in outburst together with other molecules like 
OH (hydroxyl), H I and H2 (atomic and molecular hydro- 
gen), andorganics like C2H2 (acetylene), HCN (hydrogen 
cyanide), and CO2 (carbon dioxide), giving the opportu- 
nity to investigate the conditions of the emitting gas and 
the availability and survival of these important molecules 
during the unstable early phases in disk evolution. 

2. OBSERVATIONS 

2.1. Spitzer IRS Archival Spectra 

We collected all archival spectra of EX Lupi obtained 
with the Infrared Spectrogra ph (IRS) on the Spitzer 



satellite (Houck et al. 2004). High-resolution spec- 
troscopy with the IKS is available in two different mod- 
ules with a resolving power R « 600: the short-high (SH) 
module in the wavelength range 9.9-19.6 fim and the 
long-high (LH) module between 18.7 and 37 nm. The 



IRS low-resolution modules instead cover the range 5.2- 
38 //m with a resolving power R « 60-120. The IRS 
observations of EX Lupi span the epoch surrounding its 
latest outburst in 2008. In a quiescent phase before the 
outburst two different spectra were acquired: in 2004 (PI: 
N. Evans, hereafter E04) and in 2005 (PI: G. Stringfel- 
low, hereafter S05), with both low- and high-resolution 
modules. The outburst phase was then observed in April 
2008 (PI: P. Abraham, hereafter AOS) and May 2008 (PI: 
J. Carr, hereafter C08) with the high-resolution mod- 
ules, while post-outburst spectra during a new quiescent 
phase were taken in October 2008 and April 2009 (PI: 
P. Abraham, hereafter A09) only with low resolution. In 
Table [T] all the Spitzer IRS observations of EX Lupi are 
listed with their relevant observational parameters. Since 
the resolution of post-outburst spectra (A09) is too low 
to detect the large majority of the gas emission lines 
studied in this work, we restricted our investigation to 
the high-resolution spectra taken before and during out- 
burst. More precisely, we decided to use in the present 
analysis only the spectra closer in time to the onset of the 
2008 outburst, to study the changes in emission that can 
be related to it: in the quiescent phase the S05 spectrum 
and in outburst the A08 spectrum. The investigation of 
variations/correlations in emission over a broader range 
of timescales during a quiescent period (comparing E04 
and S05) and short-term changes in outburst (comparing 
A08 and C08) may be addressed in a follow-up paper. 

2.2. Data Reduction 

We re-reduced the IRS spectra starting from the 
Spitzer Science Center S18.7.0 products and applying the 
latest version of the pipeline provided by the "Cores to 
Disks" Spitzer legacy program, described in Lahuis et 
al. (2006). This pipeline was developed to reduce IRS 
pointed observations and reach high sensitivity partic- 
ularly when observing faint objects. The extraction of 
one-dimensional spectra from the images is done using an 
optimal extraction, where the IRS point spread function, 
defined using sky-corrected high signal-to-noise (S/N) 
calibrators, is fitted on good pixels only (i.e. excluding 
known bad/hot pixels). Dither positions are first com- 
bined and then the spectrum is extracted, which gives 
the best resultant S/N. Moreover, this method provides 
an estimate of the local sky contribution, which was sep- 
arately observed with Spitzer only in the case of AOS. 
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Fig. 1. — Comparison between the quiescent (S05) and the outburst (AOS) phases in the Spitzer spectra of EX Lupi: the SH module is 
shown in the top and the LH module in the bottom. The quiescent spectrum is shifted for comparison in both plots. OH transitions are 
marked with dashed lines, while other species (C2H2, HCN, CO2, H I, H2, [Ne 11]) are marked with dotted lines. All unmarked lines are 
identified as water rotational transitions. 
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The EX Lupi spectra, like all IRS spectra of bright 
point sources, are affected by systematics that reduce 
the maximum S/N achievable with long exposures (e.g. 
uncertainties in pointing correction, flux-dependent cali- 
bration, etc.). These uncertainties are difficult to account 
for and are not included in the errors on individual pixels 
propagated through the reduction pipeline. In our anal- 
ysis we chose to derive our own estimate of the noise as 
follows. We measure in each spectrum the dispersion of 
pixels from a baseline (first-order polynomial) fitted to 
the continuum. Unfortunately, given the spectral resolu- 
tion of Spitzer, a clean continuum is difficult to find and 
often we can constrain only a pseudo-continuum given by 
weak emission lines. We use the HITRAN database to 
avoid ranges where the emission fr om o bvious molecules 
is likely to be strong (see Section 3.1 1. The pixel dis- 



persion with reference to the baseline is then chec ked for 



gauss ianity with the Shapiro- Wilk test (Shapiro & Wilk 



19651. We consider four spectral ranges in SH (11.95- 
n:^15.44-15.6, 17.45-17.7 and 19.05-19.2 fj.m) and four 
in LH (24.15-24.55, 25.45-25.85, 30-30.2 and 31.3-31.65 
fim) as an attempt to account for the variations of noise 
with wavelength. We discard the rms estimated in a 
given range if the probability from the Shapiro- Wilk test 
that the data are consistent with having been drawn from 
a Gaussian is less than 10%. As a result, in quiescence 
(S05) we find an rms of ~3 and ^9 mJy in SH and LH 
respectively, while in outburst (AOS) we find ~8 and ^-^20 
mJy respectively. We assume these estimates as the flux 
uncertainty on individual pixels in the two IRS mod- 
ules separately, and we use them i n the derivation of line 
flux errors as described in Section 13.21 On the one hand 
these should be regarded as conservative estimates as the 
molecular emission in the chosen ranges is expected to be 
small, but weak features (as well as chemical species that 
are not considered) may still contribute to the pixel dis- 
persion and thus artificially increase the noise. On the 
other hand, if high systematics strongly affects Spitzer 
spectra our estimates might still be und e restimating the 



true uncertainty (cf. Salyk et al.| 2011( |Carr & Najita 



2011 1. However, we have some mdications that our esti- 



mates of the noise are con servative rather than optimistic 
(see Sections 3.2 and 5.1). 



3. SPECTRAL LINE ANALYSIS 

3.1. Molecular and Atomic Emission 

The identification of molecular emission features in the 
Spitzer observations of EX Lupi is based on the compar- 
ison between the observed spectra and synthe tic models 
generated using the HITRAN 2008 database (Rothman 
et al.|2009 ). The MIR spectra of EX Lupi are dominated 
by a forest of water emission line^ placing this source 
within the sampl e of "wet" T Tauri disks recently re 



vealed by Spit zer ( Carr 
[2008, Pontoppidaii et al. 



Najit a|2008ll2011[ [Salyk et al 
, 2010a) . H2O lines are detected 



in both the quiescent and the outburst phases of EX Lupi 
and their measured flux with respect to the continuum 
increases in outburst over the entire spectral range cov- 
ered by Spitzer, on top of the continuum level that is 

® In the text we distinguish "lines" from "transitions" . The 
former refers to the unresolved emission features observed in the 
spectra, while the latter to the individual molecular transitions 
that produce the unresolved observed features. 




1.6 18.7 18.8 18.9 
wavelength {fJ-m) 



28.8 28.9 29.0 29.1 
wavelength {fJ-m) 



Fig. 2. — Cross-ladder OH transitions detected in EX Lupi in 
outburst. Each plot shows, from top to bottom, the outburst AOS 
spectrum, the quiescent SOS spectrum (shifted in flux for compar- 
ison), and an LTE model for OH with Tex ~ 600 K and N^^l ~ 3 
xlO^^ cm~^. In the OH model, which includes both intra-ladder 
and cross-ladder transitions, the position of cross-ladder transitions 
is indicated with vertical dashed lines. In each plot, observed emis- 
sion lines that do not match the OH model shown in the bottom are 
attributed to H2O, apart from the 14.97 ^m feature in quiescence 
due to CO2. 

^4 times higher than in quiescence (see Figure [T]) . The 
water vapor emission is composed of unresolved pure ro- 
tational transitions in the ground vibrational state with 
upper level energies (Eu) of a few 1000 K. The com- 
plex excitation structure of water makes its observational 
study a real challenge, especially at low spectral resolu- 
tion. Even small spectral ranges in the infrared can be 
finely populated by rotational transitions from different 
Eu, which are resolved only with a spectral resolution 
much hig her than that provided by S pitzer (see e.g. Fig- 



ure 5 in Pontoppidan et al. 2010a). The criticality of 



such a limitation on the derivation of the emitting gas 
properties will be addressed in Sections [4| and [5] 

In addition to water, the MIR spectra of EX Lupi re- 
veal a variety of other molecules and atoms with a no- 
ticeable difference between the two activity phases of the 
star. In quiescence the rovibrational Q-branches of C2H2 
at ^-^13. 7 /zm, HON at ^-^14 /im, and CO2 at ~15 are 
clearly detected, while a weak contribution from OH is 
identified. In outburst organic molecules are no longer 
seen, while strong H i 7-6, H i 14-9, H2 S(2), and sev- 
eral new OH lines appear and increase the confusion with 
water emission, especially in the SH module. The OH 
emission consists of rotational transitions in the ground 
vibrational state with E^ from ^600 up to ^20,000 K 
and is detected in outburst throughout the entire spec- 
trum. These transitions are mainly intra-ladder transi- 
tions belonging either to the ^H3/2 or the ^Hi/2 Pure 
rotational ladders, but several cross-ladders transitions 
are also identified, indicating a remarkable population 
of intermediate and low E^ (see Figure [2] and Table [2]). 
To our knowledge, this is the first time that cross-ladder 
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TABLE 2 
OH Detections in EX Lupi. 



A (/jm) 


Transitions 


Eu (K) 


A^i (s-i) 


S05 


AOS 


11.06 


1/2 (57/2 -> 55/2), 3/2 (59/2 -5> 57/2) 


21,300 


1300 


- 




12.65 


1/2 (47/2 45/2), 3/2 (49/2 47/2) 


15,100 


850 


V 


- 


13.07 


1/2 (45/2 43/2), 3/2 (47/2 45/2) 


13,900 


780 


- 


V 


13.55 


1/2 (4o/z — ^ 4i/Zj, v)/Z — ^ 4o/2) 


12,800 


690 




/ 


14.07 


1/2 (41/2 -s> 39/2), 3/2 (43/2 -s> 41/2) 


11,800 


620 






14.65 


1/2 (39/2 ^ 37/2), 3/2 (41/2 -s> 39/2) 


10,800 


540 


V 


if 


15.00 


1/2 (17/2) 3/2 (15/2) 


2400 


0.2 






16.66 


1/2 (15/2) -s> 3/2 (13/2) 


2000 


0.2 






18.73 


1/2 (13/2) 3/2 (11/2) 


1550 


0.2 




if 


18.85 


1/2 (29/2 -> 27/2) 


6250 


250 






23.10 


3/2 (25/2 -s> 23/2) 


4100 


130 






23.22 


1/2 (23/2 -s. 21/2) 


4150 


130 


? 




27.42 


3/2 (21/2 19/2) 


2900 


80 




i 


27.67 


1/2 (19/2 17/2) 


2960 


80 


? 




28.94 


1/2 (7/2) -5> 3/2 (5/2) 


620 


0.1 


? 





Note. — At the spectral resolution of the Spitzer IRS, all observed OH lines are blends of doublets of intra-ladder or 
cross-ladder transitions belonging to the two pure rotational ladders ^113/2 and ^Tl\/2 in the ground vibrational state. In this 
table we provide in parentheses the upper and lower J-level of each doublet. Cross-ladder doublets are reported as a single 
feature e ven when two lines can be distinguished (see Figure [2|, and in the rotation diagram analysis we use the unified flux 
(Section 5.2 1. The wavelengths we report are approximateoat the center of each feature. Molecular data are taken from 
the HITKAN database. The A^i reported here are the sums over the individual transitions blended in each observed line. 
The last two columns on the right indicate the detection in the quiescent (S05) and outburst (AOS) spectra (S/N > 2a). A 
question mark shows lines that are not unambiguously identified as OH. 



transitions are ident ified in disks, espec ially those with 
Eu ^ 1000-2000 K. |Tappe et~aL| ( |2ci08[ ), in the HH211 
outflow, were the first to ever report the detection of OH 
transitions from levels as high as ^28,200 K, and they 
also identified two lower- cross-ladder transitions (one 
of which, at 28.94 ^m, is also detected by us in the EX 
Lupi outburst spectrum). In EX Lupi, a dditional OH 
lines f rom as high as those found by |Tappe et al 



in quiescence is based on two lines in the SH module, 
while in the LH module it is highly dubious. 

Other chemical species that have been identified in the 
MIR in other T Tauri stars may additionally contribute 



(2008) might be present in the outburst spectrum, but 



their detection is hindered by the strong silicate feature 
at ^10 fim. We will come back to the analysis of OH 
emission, and the importance of the simultaneous detec- 
tion of low- and high-i^u lines in outburst, in Section 

In Spitzer spectra all the chemical species mentioned 
above are blended with H2O lines. This, in addition to 
the above mentioned blending of transitions from dif- 
ferent energy levels, hinders not only the molecular line 
identification but also the derivation of gas properties. 
Excluding the organic molecules, whose Q-branches are 
broader than the gas-phase water lines, it can be mislead- 
ing to identify molecules from only 1-2 emission lines ob- 
served where their emission is expected to be, unless such 
lines prove to be clearly separated by other species. The 
case of OH is particularly good as an example. In out- 
burst we can identify all the expected lines from the two 
rotational ladders in the ground level, thus the molecule 
is certainly detected. In quiescence, instead, only a few 
lines are unambiguously identified, at 12.65 and 14.65 
/im, which are the most spectrally separated by H2O 
lines. All other tentative OH detections are confused 
with other molecules. The case of OH lines at longer 
wavelengths (the LH module) is particularly difficult to 
judge, especially beca use of the unconstrained H2O emis- 
sion (see Section 5.1 
able models for the! 
in LH in the S05 spectrum can be explained as H2O 
emission alone. Therefore the (partial) detection of OH 



to the em ission observed in EX Lupi (see e.g. Najita et 
al. 2010). However, in the present paper we focus our 
analysis on H2O and OH, which largely dominate the 
MIR emission that we observe in EX Lupi. We assume 
that the MIR emission from these two molecules origi- 
nates in the innermost disk atmosphere (~1-10 AU from 
the central star) , as proposed fo r othe r systems in similar 
recent studies (Carr & Najita 2008, 2011; PontoppidanJ 
et al.pOlOa; Salyk et al. 2011) and confirm ed by higher- 
resolution observations toi a fe w objects (Pontoppidan 
eraL]|2010b| [Fedele et al?][20Tl] ). 



, However, consideration of rcason- 
2O emission suggests that the data 



3.2. Estimation of Emission Line Fluxes and Errors 

Given the need for reliable tracers to probe the emis- 
sion variability, we chose to use only the most trustable 
emission lines and to measure their flux and error us- 
ing a robust method based on our rms estimate. We 
restrict the analysis to strong and/or isolated emission 
lines which we believe are the least affected by artifacts. 
For instance, the lines in spectral ranges where different 
spectral orders overlap, which are affected by drop of sig- 
nal, and those where the continuum or pseudo-continuum 
is more uncertain are excluded from the analysis. In Ta- 
ble [3] we list the emission lines considered in this paper, 
indicating the spectral range of each observed feature 
(including the nearby continuum), the chemical species 
identified as contributing to the emission, and the flux 
measured in quiescence and outburst. 

Our method partly resembles techniques previously 
utilized for the measure of hne fluxes in Spitzer spec - 



tra (e.g. Pontoppidan et al. 2010a Najita et al. 2010t 



an d for the estimate of flux uncertamties (Pascucci et 
allpOS). We developed our method in IDL building on 
publicly available routines as follows. For each observed 
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TABLE 3 

Flux Estimates of Emission Lines in the IRS Spitzer Spectra of EX Lupi. 



*^T^£i/^"t" T'?^ 1 R Ti rrj3 III m i 


T ,1 Tl r II TVl 1 




(S05) 


Outburst (A08) 


11.00-11.10 


11.06 






OH 


2.63 ± 0.51 


12.22-12.50 


12.28 


OH+H2O 


1.22 ± 0.55 


H2+OH+H2O 


2.87 ± 0.43 


12.22-12.50 


12.37 


H2O 


1.51 ± 0.46 


H 1+H2O 


6.78 ± 0.59 


12.46-12.72 


12.51 


H2O 


1.23 ± 0.25 


H2O 




12.46-12.72 


12.58 


H2O 


0.62 ± 0.23 


H 1+H2O 


3.71 it 0.61 


12.46-12.72 


12.64 


OH 


0.82 ± 0.26 


OH 




12.67-12.90 


12.81 


[Ne II ]+OH 


0.66 ± 0.28 


[Ne II ]+OH 


2.59 ± 1.13 


12.95-13.15 


13.06 


H2O+OH 


0.73 ± 0.38 


OH 


1.33 ± 0.51 


13.40-13.60 


13.43 


H2O 


0.90 ± 0.14 


H2O 


3.27 ± 0.53 


13.40-13.60 


13.49 


H2O 


0.39 ± 0.16 


H2O 


0.86 ± 0.31 


13.40-13.60 


13.55 


OH 




OH 


1.51 ± 0.44 


13.50-13.80 


13.70 


C2H2+H2O 


2.76 ± 0.29 






13.75-14.12 


14.00 


HCN+H2O 


7.48 ± 0.63 






13.95-14.30 


14.06 


OH 




OH 


1.22 ± 0.37 


13.95-14.30 


14.20 


H2O 


1.38 ± 0.30 


H2O 


1.71 ± 0.56 


14.43-14.61 


14.51 


H2O 


1.87 ± 0.20 


H2O 


2.08 ± 0.46 


14.55-14.85 


14.65 


OH 


1.10 ± 0.39 


OH 


3.51 ± 1.19 


14.80-15.05 


14.95 


CO2+H2O 


3.20 ± 0.29 


H2O 


3.36 ± 0.77 


14.93-15.10 


15.00 






OH 


1.35 ± 0.66 


15.67-15.87 


15.74 


H2O 


1.48 ± 0.17 


H2O 


2.04 ± 0.45 


15.85-16.10 


16.00 


H2O 


1.43 ± 0.22 


H2O+OH 


4.34 ± 0.95 


16.05-16.18 


16.11 


H2O 


0.85 ± 0.19 


H2O 




16.53-16.80 


16.66 


H2O 


1.47 ± 0.18 


H2O+OH 


2.87 ± 0.76 


16.73-16.99 


16.89 


H2O 


1.34 ± 0.22 


H2O+OH 


5.76 ± 0.84 


17.03-17.19 


17.11 


H2O 


1.01 ± 0.14 


H2O 


1.61 ± 0.66 


17.14-17.31 


17.22 


H2O 


1.26 ± 0.15 


H2O 


3.36 ± 0.60 


17.28-17.46 


17.36 


H2O 


2.14 ± 0.28 


H2O 


2.70 ± 0.39 


18.65-18.80 


18.73 






OH 


1.45 ± 0.55 


18.73-19.11 


18.85 






OH 


2.17 ± 0.61 


18.73-19.11 


19.00 


H2O 


1.24 ± 0.11 


H2O 


3.37 ± 0.97 


20.68-20.90 


20.79 


H2O 


2.85 ± 0.28 


H2O 


2.78 ± 0.74 


21.70-22.04 


21.84 


H2O 


5.51 ± 0.44 


H2O 


6.52 ± 0.57 


22.82-23.37 


23.10 






OH 


7.26 ± 2.87 


22.82-23.37 


23.22 


H2O 


0.50 ± 0.19 


OH 


5.64 ± 1.64 


27.20-27.88 


27.42 






OH 


6.43 ± 1.43 


27.20-27.88 


27.67 


H2O 


1.21 ± 0.35 


OH 


4.58 ± 0.82 


28.81-29.07 


28.94 


H2O 


1.23 ± 0.22 


OH 


1.45 ± 0.56 


29.72-30.10 


29.84 


H2O 


1.66 ± 0.27 


H2O 


9.39 ± 1.19 


30.40-31.05 


30.49 


H2O 


1.05 ± 0.26 


H2O 


6.88 ± 0.92 


30.40-31.05 


30.88 


H2O 


3.12 ± 0.25 


H2O 


4.61 ± 0.51 


32.78-33.40 


33.00 


H2O 


5.03 ± 0.36 


H2O 


8.62 ± 0.82 



Note. — We give t he sp ectral range of the observed emission features where we apply our method to derive line fluxes 
and errors (see Section 3.2 1. For each range we indicate the identified species most contributing to the measured flux. In 
the cases where it was possible to distinguish different lines within a feature (e.g. see Figure |3|, we report the components 
separately. The line flux is given in units of 10"^'' erg s~^ cm~'^. All fluxes reported in the table have S/N > 2a. Where an 
entry in the table is omitted, emission lines are not clearly identified and/or the measured flux is consistent with zero. 



emission line we select a spectral range including at least 
three pixels on each side for the continuum (typically 
a few tenths of a micron). Where a continuum cannot 
be found in between two or more lines because of their 
vicinity, we fit them together. We fit the data using 
a least-squares Levenberg-Marquardt algorithm (Mark- 
[ward t 2009 ) with a baseline plus one Gaussian function 
for each line considered, with width and area as free pa- 
rameters. Since we fit the continuum locally, we always 
assume a first-order polynomial for the baseline even 
though a higher order could in some cases slightly im- 
prove the fit. Each line flux is computed as the area of 
the locally continuum-subtracted best-fit Gaussian. The 
observed emission line(s) must be well represented by one 
(or more) Gaussian functions, and we set an acceptance 
threshold to the goodness of the fit to Xred ^ In 
most cases the Xred niuch lower than that limit, sup- 
porting that our noise estimate is indeed conservative 
(see e.g. Figure [S]). This method provides the possibil- 
ity to measure the individual line fluxes even in the case 



of complexes composed of several blended lines from the 
same or different chemical species (an example is shown 
in Figure [3|, unless the different lines are overlapping too 
much to be distinguished, as it is often the case in the LH 
module. This approach is not used for the broad emission 
features from C2H2, HCN, and CO2, because the shape 
of the Q-branch, which is at least partially resolved by 
the IRS, is not Gaussian. We therefore fit only a base- 
line continuum to nearby regions as described above and 
then integrate the flux below the continuum-subtracted 
data. 

The 1-(T error on each line flux is then estimated using 
the following Monte-Carlo based approach. In each of 
the spectral regions considered above, we add normally 
distributed noise to the observed flux of individual 
pixels, that we take as the mean. As the noise to be 
added we use the rms calc ulated from the supposed 
line- free regions (Section |2.2[ ). The procedure is repeated 
1000 times and we measure the line fluxes from each 
realization exactly as described above. The error on the 
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Wavelength (//m) 

line 2 



Flu< {erg s"' cm-=) 




Fig. 3. — Measure of line flux and error for emission lines in the 
SOS spectrum. In the first plot (top, left) the best-fit curve (first- 
order polynomial -|- three Gaussians) is overplotted to the original 
data, where the error on in divid ual pixels is set to the measured 
rms dispersion (see Section 2.2 i. The Xred °^ best fit is 0.2. 
Such a low value suggests tnac our rms estimate is conservative. 
The error on each line flux is set to a robust standard deviation of 
the distribution of 1000 repetitions of the flux measure (see Section 
|3.2|l. The first two lines are H2O lines, while the third is OH. 



line flux is then set to a robust estimate of the standard 
deviation of the distribution of its 1000 realizations 
(see Figure [3]). With this method we take into account 
the noise on single pixels (given by our estimate of the 
rms) as well as the uncertainty on the local fit to the 
continuum. 



4. CONSTRAINING GAS PROPERTIES 

4.1. An LTE Approach 

A simple approximation utilized in observational as- 
trochemistry is to assume the emitting gas to be in local 
thermal equilibrium. In LTE the population of the ex- 
cited levels simply follows the Boltzmann distribution, 
i.e. all levels are thermalized to the kinetic tempera- 
ture (Tk) of the gas. A transition between two levels can 
be described as to be in LTE when collisions dominate 
the excitation/de-excitation, i.e. when the density of the 
ambient gas is higher than the critical density of the 
transition (jidit = Aui/Cui, where A^i is the Einstein- 
A coefficient and C„/ the collision rate). If all transi- 
tions are in LTE, the kinetic temperature and the col- 
umn density of the gas can be derived from the observed 
flux of spectrally-resolved transitions. The techniques 
that are most frequently used to derive gas properties 
are based on the rotation diagram analysis (e.g. ,Gold-, 
smith & Langer 1999 herea fter GL99) or the dir ect fit 

here- 



ot transition intensities (e.g. Meijerink et al. 



2009 



after ME09). In the extremely diverse conditions found 
in astrophysical environments, however, it is often the 
case that the ambient density is insufficient to thermalize 
some (or all) levels and the LTE description falls short 
of deriving the real gas properties. ME09 showed how 
the LTE approximation is inappropriate for circumstel- 
lar disks, as the wide range of thermodynamic condi- 
tions found in the inner ~ 10 AU of disks (gas tempera- 
tures ranging from ~ 100 to a few 1000 K, and densities 
from ~ 10'^ to 10^^ cm~'^) results in emission produced 



by a mixture of thermalized and non-thermalized levels. 
GL99 showed how the rotation diagram analysis can help 
in distinguishing the effects of optical depth and non- 
thermal excitation on the typical straight-line behavior 
of optically thin lines in LTE. However, for this type of 
analysis spectrally resolved observations are required, to 
separate different energy levels. With the spectral res- 
olution provided by Spitzer this is possible in the case 
of OH, but not water. The implementation of such a 
techniq ue in the case of OH emission will be presented in 
Section [K2I Now we focus on the method used to derive 
constraints on the excitation properties of the unresolved 
and blended water emission. 

Non-LTE models that perform line radiative transfer 
and account for the geometrical and physical struc- 
ture of disks have been e xplored to some extent (e.g. 



Pavlyuchenkov et al. 2007), but they are not yet viable 
alternatives tor constraining water MIR line emission be- 
cause of the incompletenes s and uncertainty of collision 
data (see van der Tak|20H ) . This is why the simple LTE 
approximation is still the first approach utilized when 
mo dellin g the molecula,r emission observed in Spitzer 
spectra |Carr fc Najita||2008l [20lT| |Salyk et al.] [2011). 
We therefore developed our own simple single-slab Li'E 
model to constrain the molecular emission in EX Lupi. 
The model is explained in detail in Appendix 1X1 while 
here we provide only a short description. The free 
parameters of the model are the excitation temperature 
Tox, the column density N^oi, and the emitting area A. 
In the fit results we will report the value of the radius, 
Tp (projected), of the circular area that is equivalent 
to A. The excitation and relative strength of different 
transitions are sensitive to Tex and iVmob while A only 
has the effect of a scaling factor common to all line 
fluxes. Observed and model spectra are compared using 
the individual fluxes of blended emission lines (given 
in Table pi), and a test is performed minimizing 
over the oifference between observed and model line 
fluxes as an ensemble. Despite the fact that the model 
parameters are sensitive to different properties and 
that we account for line optical depth effects, the fit 
results are degenerate in Tex, A^moi, and A. In general, 
for a given emitting area the fit is degenerate with 
higher Tox for lo wer iVmoi (and vice versa, see e.g. 
results in Section [5.1[) . In the extreme optically thick 
regime (t ^ 1) the line flux depends on the product 
A X Tox, thus only those can be derived while the 
model is insensitive to changes of Nnmi- In the extreme 
optically thin case (r ^ 1), instead, the line flux is 
proportional to A x N^^i x Tex, and we can reduce again 
to two parameters: the temperature and the number 
of emitting molecules = A x Ny^oi, yet still suffering of 
the degeneracy between Tqx and iVmoi- Ideally, three 
emission lines that are differently sensitive to variations 
in Tox and Nnmi should be enough to constrain all 
model parameters together. Water MIR transitions 
do cover a wide range of opacities and therefore an 
appropriate selection of them could serve the scope. 
However, in practice their different opacity cannot 
be retrieved unambiguously from the blended lines 
they produce in the Spitzer spectra and the degener- 
acy between model parameters is somewhat unavoidable. 
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4.2. Fitting Spitzer Spectra of H2O Emission with LTE 

Models 

Fitting the MIR unresolved water emission in EX 
Lupi using LTE single-slab models is challenging. One 
must keep in mind that such simplified models do 
not account for any disk geometry, structure, nor 
tempera ture/column-densit y profiles, which are complex 
(see e.g. iBergin et al. 2007) and in the case of EX Lupi 
might also change trom quiescence to outburst. In addi- 
tion to that, fiux contamination from different molecules 
can also affect the fits, unless different contributions can 
be individually constrained and removed (which usually 
require many "clean" lines of the same species and a re- 
liable modelling). In fact, we found that in EX Lupi 
it is impossible to get a good fit over the wide spectral 
range covered by the Spitzer IRS with a single set of 
Tex and A^moi- Yet it is possible to obtain a reasonable 
fit for sub-sets of lines, that likely share similar exci- 
tation conditions. Even in the literature there is still 
no common agreement on the most appropriate way to 
fit single-slab LTE models to th e water MIR emis sion 
from disks observed with Spitzer. Salyk et al. (2011) try 
to derive average gas properties for 48 1' 'I'auri systems 
fitting peak-to-continuum values of many emission lines 
over the entire 10-35 fiin range, but they also note that it 
is impossible to reproduce entirely the observations using 
a single-slab LTE model and that the outcome varies if 
different spectral ranges an d/or fitting methods are used. 
Carr fc Najita (2008 2011 1 instead focus on fitting lines 
observed in a small spectral range limited between 12 and 
16 fim. Applying this method to a sample of 11 T Tauri 
systems Carr & Najita (2011 1 report that while they can 
reach a very good match over this range, the same model 
overpredicts lines at longer wavelengths (16-19 /xm). Dif- 
ferent methods provide different results, so what is the 
most appropriate way to proceed and how shall we in- 
terpret different gas properties obtained fitting different 
spectral ranges? 

A viable solution comes from theoretical works. ME09 
showed that LTE models generally overpredict the line 
flux if H2O is not in LTE, as emission lines excited in 
non-LTE regimes are sub- thermally populated (see also 
[Pavlyuchenkov et al.|2007 who explored the case of other 
molecules). In their Figure 8, ME09 show that the dif- 
ference in line flux of a water transition in the ground 
vibrational level excited in LTE and in non-LTE condi- 
tions is minimal for small values of A^i- If the excita- 
tion state of the emission is unknown, a sufficient spec- 
tral resolution may allow a selection of lines to fit with 
LTE models based on the value of Aui- However, they 
also show that such transitions are widely distributed in 
wavelength (Figure 7 in ME09). Therefore, when the 
resolution is not enough to separate different transitions 
(as in the case of the Spitzer IRS) it is hard to say in ad- 
vance how well a blended line can be described assuming 
LTE excitation. If the emission is originated in a mixture 
of LTE and non-LTE conditions, which is a reasonable 
assumption for T Tauri disks (ME09), it is therefore not 
surprising that LTE models may not be able to reproduce 
the relative strength of the observed lines. Moreover, as 
ME09 remark, even finding a good fit {Xred ~ 1) a. 
handful of blended water lines does not necessarily in- 
dicate that the LTE model provides a good estimate of 



the real conditions of the emitting gas. It is clear that 
some criteria for the fitting method are needed in order 
to use LTE models and give reliable results that do not 
misinterpret the data. 

A sensible choice is to use LTE models to fit only lines 
at the shortest wavelengths, where the overprediction 
of the line flux assuming LTE is minimal (Figure 7 in 
M E09). This is c o nsistent wit h the approach utilized 



by Carr & Najita 
should check the va 



2008 



ue 



H 

of 



2011|. In this approach one 
tlu and Aui of the transitions 



blended in the observed lines, because transitions with 
high E^ii and A^i can be relevant even at short wave- 
lengths (ME09). Another possibility is to consider the 
critical density and select lines based on the lowest val- 
ues of ricrit, which are the first to be thermalized if the 
ambient gas density is high enough. However, collision 
rates of water transitions are a ccurate only with in a fac- 
tor of 3-10 depending on E^ ( Faure fc Josselin,,2008) , 
which in turn affects rtcrit and its use as a selective tool. 

4.3. Fitting Method Adopted in This Work for H2 
Emission 

In this work we chose to fit those H2O lines that are 
the least likely to be affected by non-L TE e ffects. Based 
on the arguments presented in Section |4.2[ we start fit- 
ting over the short wavelengths (12-15 /im) and then we 
explore changes in the outcomes adding/removing lines 
in the range 12-19 /im. Each time a best-fit model is 
found for a sub-set of lines, we check the E^, A^i, and 
ncrit of the transitions that contribute for at least 10% 
to th e line fluxes. For ricriti we take the collision rates 
fromlFaure & Josselin ( |2008 ) and we consider collisions 
with H2 only. We try to avoid lines that include tran- 
sitions with Eu ^ 4300 K or A^i > 1 s""'^, as in such 
cases non-LTE effects are more likely to be strongly rel- 
evant (ME09). To check our best-fit models we consider 
two factors: the value of the xled (that we require to be 
less than 1.5) and the plot of the ratios of observed and 
model line fluxes against the wavelength. More precisely, 
we expect the latter to resemble the behavior reported in 
Figure 7 in ME09 if the fitted lines are in fact the closest 
to a thermal excitatiorj^ In other words, an LTE model 
should be able to reproduce the flux of lines given by 
transitions with the lowest values of and ricritj 

while in the case of transitions with higher values of 
Aui, and Ud-it the line flux should be strictly overpre- 
dicted. 

It is clear that our approach should be considered ex- 
ploratory. Moreover, the fact that a single-slab LTE 
model fails to reproduce the observations may imply that 
the observed emission is indeed excited in non-LTE con- 
ditions and/ or that it may be produced in several disk re- 
gions with different Tex, A^moi, and emitting areas rather 
than in a single homogeneous layer. Given the limita- 
tions of using a single-temperature LTE model, we are 
able to constrain only the emission that originates in a 
portion of the disk atmosphere where the local gas den- 
sity is high enough for thermalization, and which is nar- 
row enough to be approximated with one temperature. 
However, thanks to the test over the line ratios as ex- 
plained in the previous paragraph, we are able to obtain 

^ Note that we plot the ratio of observed/LTE fluxes while ME09 
plot non-LTE/LTE fluxes. 
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Fig. 4. — Quiescence: ratios of the observed and LTE model water line fluxes for the model fits to lines in SH (left) and over the entire 
IRS spectrum (right). The model parameters are reported in the box. 
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Fig. 5. — Same as Figure [4] but showing the best-flt found using only the lines indicated by filled circles, which are the least likely to 
be affected by non-LTE effects. All other lines, shown as empty circles, are consistent with being closer to no n- LTE conditions and are 
overpredicted by the LTE model, while in LH the best-fit found in SH is clearly not appropriate. 



indications on whether the emission is mainly to be at- 
tributed to non-LTE excitation rather than to different 
tem pera tures from different regions in the disk (see Sec- 
tion 5.1 ). 



5. RESULTS OF MODELLING MOLECULAR EMISSION 

5.1. H2O Rotational Lines 

For the sake of comparison with the results obtained 
following our "LTE-fitting" method explained above in 
Section |4.3[ we first report the results obtained fitting 
lines regardless of the values of , A^i , and ricrit ■ In qui- 
escence, the results obtained from fitting water lines over 
the SH module first, and then over the entire SH+LH 
coverage, are shown in Figure H] The temperature varies 
from ~ 520 to ~ 330 K in tne two cases respectively, 
while the column density stays close to ^ 5 x 10^* cm~^. 
The probability that the data are well described by the 
models is lower than 0.001% in both cases, i.e. the LTE 
model is strongly inappropriate as the scatter in Figure 
|4] suggests. We therefore proceeded, as described above, 
with fitting sub-sets of lines at short wavelengths, and we 
found a best fit over the range 12-17 /im with T^x = 580 
K and iVmoi = 3.2 x 10^^ cm^^ (ggg Figure [5]). The xled 
is much lower than 1, suggesting that our error-bars may 
overestimate the real errors by a factor of ~ 5 or more. 



The lines used in the best fit are predominantly produced 
by transitions with Eu < 4300 K and A^i < 1 s~^, and 
their critical densities range between 10^ and a few 10^ 
cm^^ (see Table [4]). Therefore the chosen lines are the 
first to thermalize, compared to other lines in the SH 
module that instead include transitions with E^^ up to 
7000 K, Aui up to 10^ s"\ and ricnt up to 10" cm"^. 
The behavior of the line flux ratios plotted in Figure 
[5] supports that we are indeed probing the most ther- 
malized lines: while the LTE model can reproduce very 
well the fitted lines, other lines are instead strictly over- 
predicted. The emission is close enough to an optically 
thick regime (r ^ 1-10) that the emitting area is well 
constrained and the best fit is found for Vp = 0.58 AU. 
The confidence limits derived for the best-fit parameters 
are shown in Figure [8] and fisted in Table [5] 

The same procedure is applied to the outburst spec- 
trum, and the results are shown in Figures [6[ and [Tj Here 
we have an additional limitation given by the increased 
confusion of water lines with other molecules/atoms, 
mainly OH and H i at short wavelengths (see Figure [T|. 
Therefore, from the sample of lines that were used m 
quiescence we exclude those which are the most contam- 
inated by other species (see Table [3]) . Fitting over the 
SH module gives a temperature of ~~370 K and a column 



10 






T= 373 K 


N = 


.le + 20 cm"^ 




r= 0.87 AU 



T= 293 K 
N= 4.8e+18 cm"^ 
r= 2.61 AU 



wavelength (^um) 



20 25 
wavelength (^m) 



Fig. 6. — Outburst: ratios of the observed and LTE model water line fluxes for the model fits to lines in SH (left) and over the entire 
IRS spectrum (right). The model parameters are reported in the box. Water lines that are heavily contaminated by other species are not 
considered here. 
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Fig. 7. — Same as Figure [6] but showing the best fit found using only the lines indicated by filled circles, which are the least likely to 
be affected by non-LTE effects. All other lines, shown as empty circles, are consistent with being closer to no n- LTE conditions and are 
overpredicted by the LTE model, while in LH the best-fit found in SH is clearly not appropriate. The lino at ~ 13.43 fim is strongly 
underpredicted by this model and is excluded from the fit. 



density of ~ 10 cm" , while over the entire spectrum 
we find a temperature of ~ 300 K and a column density 
of ^ 5 X 10^^ cm^^. Again, the model is not appropri- 
ate to reproduce the observations (probability lower than 
0.001%), as the scatter in Figure p] shows. The best-fit 
model using instead only lines closer to thermalization 
gives a good fit with T^^ ~ 600 K and Ny^oi ~ 3 x 10^'' 
cm~^, but the emission is more optically thin than in 
outburst (r ^ 0.1-1) and it is sensitive to changes in the 
emitting area only for low values of Tp (i.e. high values 
of Nyno\, when the emission gets optically thick). The 
emission is found consistent with Vp > 0.6 AU (68% con- 
fidence), but for Tp > 2 AU it is so optically thin (t < 
0.1) that the result is completely insensitive to a change 
in emitting area. Assuming that OH and H2O emission 
probes the same disk radii (as recentl y shown from very 
high- resolution near-infrared lines by |Doppmann et al. 



2011 1, we choose to fix the emitting area for H2O in out- 



burst to Tp = 1.3 AU ( the best-fit value of low-E^ OH 



in outburst, see Section 5.2). Then we derive confidence 
limits on T^x and N^oi tor this fixed radius, which we 
take as our reference case (see Table [s] and Figure [s]) . 
The opportunity of having the same emitting area for 



OH and H2O is that we can easily compare their rela- 
tive change in terms of Tex and N^oi only. In the case 
that the assumption was found to be wrong and H2O had 
instead the same radius found in quiescence (0.58 AU), 
the model gives Tex « 700 K and A'^moi « 1 x 10^^ cm'^ 
to explain the increase in line fluxes. In our reference 
case, fixing rp to 1.3 AU, the lines considered in the best 
fit have mostly transitions with ranges of E^, A^i, and 
'T-crit consistent with the quiescent case, even though the 
best-fit model shows a relevant contribution from tran- 
sitions having higher values of Auu and n„it (see Table 
[4]) . The behaviour of line ratios discussed above is once 
again confirmecp"] (Figure [t]) . However, it is clear that 
the sub-set of lines we end up fitting in outburst is more 
affected by non-LTE excitation, and that the Tox derived 
is therefore probably only a lower limit to the real gas 
temperature of the observed emission. 

In fact, Tox always lowers when we include in the fits 
lines that may be more affected by non-LTE according 

A noticeable exception is the line at ^ 13.43 ^tm, which, 
despite being one with the lowest Eu and A^i, is underpredicted 
by the best-fit model. Contamination by another species is a viable 
explanation. 
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TABLE 4 

Parameters of HtO Lines Used to Constrain the Emission Observed in EX Lupi. 



Line (^m) 



Transitions (/im) 



Eu (K) 



ncrit (cm 3) 



Quiescence (SOS) 



12.51 


12.51 


4133 


1.50 


4(9) 


0.6 




12.52 


4133 


1.50 


3(9) 


1.7 


14.51 


14.49 


3234 


1.22 


3(9) 


2.1 




14.51 


3233 


1.21 


2 (9) 


6.2 




14.54 


2605 


0.15 


3(8) 


1.6 


15.74 


15.74 


2824 


1.09 


2 (9) 


11.4 


16.89 


16.89 


2698 


0.70 


2 (9) 


8.2 




16.90 


2125 


0.26 


5 (8) 


7.5 



Outburst (AOS) 



13.49 


13.48 


5499 


8.61 


2 (10) 


0.1 




13.50 


3341 


0.49 


1 (9) 


0.2 


14.20 


14.18 


2876 


0.30 


6(8) 


0.2 




14.21 


3951 


3.35 


8(9) 


0.5 


14.51 


14.49 


3234 


1.22 


3(9) 


0.2 




14.51 


3233 


1.21 


2 (9) 


0.6 




14.54 


2605 


0.15 


3(8) 


0.2 


15.74 


15.74 


2824 


1.09 


2 (9) 


1.1 



Note. — For each observed line, we report the transitions that, according to the best-fit model, contribute for more than 
10% to the observed flux. The value of the opacity of each transition is dependent on the best-fit model (in outburst we 
report here the model with r-p = 1.3 AU). All other parameters are instead intrinsic of each single transition. Molecular data 
are taken from the HITRAN database, apart from the collision rates which are taken from Faure & Josselin (20081. Critical 
densities must be multiplied by 10", where the exponent a is shown in parentheses. 



to the criteria explained in Sections |4.2| and |4.3[ It is 
very interesting that a similar effect is seen in comparing 
two recent studies. Despite the similarity in the samples 
of T Tauri stars considered in t he two pape r s, the aver- 
age gas temperature reported in Salyk et al. (2011 1, who 
include in their fits lines a t wavelengths li p to 3 5 /im, 
is » 25% lower than what Carr & Najita (2011) find, 
who instead limit their fits to the 12-16 /xm range. A vi- 
able explanation is that the different results are due to a 
different contribution from non-LTE excitation given by 
the different choice of the spectral ranges to fit. This is 
supported by the fact that, using the outburst spectrum, 
our result obtained from fittin g lines over the entir e spec- 



trum is consistent with what Salyk et al. (2011) found 



for EX Lupi in outburst, while when we restrict the fits 
to lines closer to LTE the results are consistent with the 



average found by Carr & Najita (2011 ). As already said, 
if transitions at long wavelengths sutter more from non- 
LTE effects (i.e. are sub-thermally populated), then the 
Tox found using an LTE model fitted to such transitions 
will be lower than if they were properly modelled in non- 
LTE. The fact that a lower temperature is always found 
when lines at long wavelengths are included strongly sug- 
gests that the MIR emission we observe is indeed from 
non-LTE excitation, and that only some lines can be ap- 
proximated as being in LTE. 

However, non-LTE excitation is not alone in affecting 
MIR molecular emission. In fact, the best fits we obtain 
from fitting lines closer to LTE are clearly not appropri- 
ate for the emission lines detected in LH: line fluxes are 
strongly undcrpredicted, contrarily to what is found in 
the SH module (see Figures [5] and [7]). Lines at wave- 
lengths longer than 20 /im are predominantly produced 
by lower Eu than at shorter wavelengths, typically ^ 
900-1800 K. Therefore an LTE solution with Tex ~ 200- 
300 K could in principle account for such low levels. How- 



ever, the strict overprediction of line fluxes that we find 
using lines in the SH module cannot be found using lines 
in LH even if we assume a low temperature. A reason 
for that can be seen in that the A^i are instead gener- 
ally higher {A^i > 10 s~^) than at short wavelengths, 
and higher is in turn the effect of non-LTE excitation. 
This is again consistent with Figure 7 in ME09: at wave- 
lengths longer than 20 /im the ratios of non-LTE and 
LTE fluxes span the whole range from 1 to 0, such that 
no line observed with the Spitzer resolution can be ap- 
proximated as in LTE. Nonetheless, non-LTE excitation 
alone is not able to account for the observed fluxes at 
long wavelengths when we assume the same temperature 
found at the short wavelengths. Given the lower E^, 
what we observe in the LH module is likely produced 
by water vapor emitted from a colder region in the disk, 
probably at larger radii. 

5.2. OH Rotational Lines 

OH rotational transitions from different levels are not 
as intermingled in wavelength as for the H2O molecule, 
and the individual lines observed with Spitzer are com- 
posed of doublets from a very narrow range of E^ and 
Aui (see Table [2| . This enabled us to perfor m a rota- 
tion d i agrarn . analysis, using the description in ILarsson] 



et al. 



( |2002[ ). In Figure I9| we show the rotation dia- 
gram based on the OH lines detected in outburst. We 
assume that their measured fluxes are entirely given by 
OH emission, and we exclude the ~ 16.66 /im doublet 
which is strongly blended with water (see Figure |2|. It 
is apparent from the rotation diagram that the detected 
lines are not well approximated with a single straight 
line, that would be expected if the emission was opti- 
cally thin and in LTE. GL99 showed how deviations can 
be due to different factors: LTE/non-LTE excitation, op- 
tically thin/thick emission, the complexity of the molec- 
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N (cm ^) 

Fig. 8. — Confidence limits on Tex and N^^-^^i estimated for H2O in quiescence (SOS, left) and in outburst (AOS, middle), and for the 
low-i?„ OH co mpon ent in outburst (right). In the case of water in outburst we show the contours for a fixed emitting area (with rp=1.3 
AU, see Section jsTT} . Contours show the 68%, 95% and 99% two-dimensional confidence regions, while the best-fit locus is marked with a 
star. The degeneracy between the two model parameters is apparent, giving higher/lower temperatures for lower/higher densities. 

2 [ 



ular structure, and a multi-temperature emission. The 
interpretation of the emission we observe in EX Lupi is 
not straightforward because we may have a combination 
of effects. First of all, the OH molecule does not belong 
to the simple linear molecules studied in GL99, given its 
ground-state double-ladder excitation structure. Second, 
we very likely observe emission that can only partially be 
approximated to LTE, as in the case of H2O. Third, the 
opacity of emission lines is not known a priori and de- 
pends on Tox and A^moi- Unambiguously distinguishing 
between these different effects is challenging. We propose 
here a viable interpretation. 

Given the much smaller transition probabilities of the 
cross-ladder compared to the intra-ladder transitions (see 
Table [2]) , the former are the most likely to be optically 
thin. A linear fit to the cross-ladder fluxes alone finds 
Trot ~ 950 ± 460 K, but we note that a reasonable fit 
could also include lines up to Eu ^ 6000 K. If we extend 
the linear fit to include such lines we find Tiot ~ 620 ± 50 
K (see Figure |9]). It is instead apparent from the dia- 
gram that such a low rotational temperature is certainly 
not adequate for higher-i?„ lines. When fitting sepa- 
rately over the high-_Ei, range {E^ > 10, 000 K) we find 
Trot ~ 9500 ± 5200 K (see Figure |9]). Two components 
can explain the observed line fluxes even when we ac- 
count for the o ptica l depth using our LTE model pre- 
sented in Section 4.1 A model fit to lines with E^ < 6000 
K finds Tex = 5g(rK and N^oi = 2.8 x 10^^ cm^^ ^^^^ 
an emitting area with Tp — 1.30 AU. The emission is 
optically thick (although the cross-ladder transitions are 
still optically thin) and well characterized (see Figure |8] 
and Table [5]) , so that we take its area to be t he same 
for water in outburst as assumed in Section [5.11 On the 
other hand, also the high-_Ei, lines {Eu > 10, 000 K) can 
be reproduced reasonably well by a single-temperature 
LTE model with Tox ~ 10, 000 K and N^oi « 1 x lO^"* 
cm~^ (fixing rp to 1.3 AU), the emission is very optically 
thin and cannot reproduce line fiuxes with Eu < 6000 
K. Despite the success of these two separate LTE so- 
lutions for different ranges of Eu (see Figure l9| we do 
not conclude that the emission is in LTE. A few tran- 
sitions lying on a straigth line can also be mimicked by 
a quasi-thermal behavior (GL99), and finding different 
temperatures when different ranges of Eu are fitted sug- 
gests non-LTE excitation. In conclusion, this analysis 
shows that the emission in EX Lupi in outburst can be 
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Fig. 9. — Rotation diagram of OH lines detected in outburst 
(see Table [2I, reported as black squares with error-bars. L inea r 
fits suggest two separate emission components (see Section |5.2| , 
and are shown with red/blue dotted lines. Two LTE models can 
reproduce the line fluxes of each component separately (red/blue 
crosses), but neither one is able to reproduce all lines together (red 
crosses fall outside the plot in the high-iJ^ range). Consideration 
of the two detections in quiescence (empty circles with error-bars) 
suggests a lower column density and a temperature consistent with 
outburst for the component given by high-Eu lines. 

explained by two LTE components, one warm (T ^ 600 
K) and one hot (T > 4000 K). It is nonetheless fikely that 
a two-temperature LTE description is simply an approx- 
imation, and that many factors play a role in producing 
the scatter/curvature that we see in the rotation dia- 
gram: optical depth at low Eu, non-thermal excitation 
at high Eu (in fact for Tgas ^ 10, 000 K thermal disso- 
ciation would become important), and/or a range of gas 
temperatures instead of only two. 

We briefly analyze also the quiescent case. The high- 
Eu component, with only two detected lines, suggests 
a comparable temperature and a lower column density 
as compared to outburst (see Figure [9). As mentioned 
above, the low-Eu component is instead not detected. If 
we assume that the flux from the two dubious detections 
at - 23.22 and 27.67 (Figure IT] and Table [2]) is to 
be attributed to OH alone and we tit these lines, we find 
a low-i?„ OH column density of ~ 0.7 x 10^^ cm~^ (as- 
suming the same Tex found in outburst, 590 K, and the 



same emitting area found for water in quiescence, 
rp = 0.58 AU). 



with 
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high-E„ OH 




Fig. 10. — Quiescent spectrum (SOS) and LTE models derived 
for H2O (blue) and OH (cyan), plotted together with the models 
assumed for HCN (green), C2H2 (purple) and CO2 (orange). The 
sum of all model contributions from different molecules is shown 
in black and shifted for comparison. A continuum is not added to 
the model. For the three organic molecules we assume the same 
emitting area found for water emission (rp = 0.58 AU), and model 
parameters consis tent with the range of values derived for other T 
Tauri systems by |Carr fc Najita| |2008| |2011| , namely Tex = 750 
K and 7Vm„i ~ 2-7 xlO" cm--=. 
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Fig. 11. — Same as Figure [To| but showing the outburst (AOS) 
spectrum. The LTE model derived for \ow-Eu OH emission is 
added in grey (two lines can be seen at ~15 fim). Models for 
organics are plotted using the same parameters assumed in quies- 
cence. It can be seen that in outburst the Q-branches of C2H2 at 
~13.7 fim, HCN at ~14 fim, and CO2 at ~15 ^m are not as strong 
as in quiescence. Remarkable emission lines from other identified 
molecules/atoms are labelled. Unidentified species may be respon- 
sible for other features, like those at ~ 13.8 fim. 



5.3. C2H2, HCN and CO2 Rovibrational Branches 

The Q-branches of C2H2, HCN and CO2 produce 
broad blended features in the Spitzer spectra, contribute 
to each other's flux, and are also contaminated by H2O 
and OH lines. The shape of the Q-branch, as it is seen 
with the Spitzer spectral resolution, depends on both the 
temperature and the column density of the emitting gas. 
Therefore, our fitting method is not sufficient to con- 
strain the emission from organics, as the line flux alone 
is very highly degenerate in Tox and iVmoi- A careful 
analysis of the unresolved emission from these molecules 
requires at least that all other chemical species are well 
constrained and subtracted from the data, and that the 
Q-branch profile is accounted for in the model fitting. 
T wo attempts in this dire ction have been recently shown 
by |Carr fc Najita| ( |20TT[ ) and jSalyk et al.| pOTl] ). The 
constraints on 7cx and N^oi they are able to derive are 
loose such that a characterization of the change in or- 
ganic emission in the case of EX Lupi between quies- 
cence and outburst would be very hard and uncertain. 
Therefore in this paper we refrain from any kind of so- 
phisticated analysis of the organic emission, and simply 
provide a rough estimate of the change observed between 
quiescence and outburst in terms of a variation in column 
density. 

In Figure [TO] we show the quiescent S05 spectrum and 
LTE models for C2H2, HCN, and CO2 assuming parame- 
ters in the ranges deriv ed for other T Tauri disks by |Carr| 
k Najita] ( |2008[ |201l| ). Namely, for EX Lupi we assume 
Tex = 750 K and TV^oi ~ 2-7 xlO^^^ cm"^ while we take 
the emitting area from water emission (rp = 0.58 AU). It 
can be seen from the figure that the emission observed in 
EX Lupi does not sho w a ny remarkable difference from 
the models. In Figure 11 it can be seen that the models 
assumed for organics in quiescence are instead not con- 
sistent with the outburst AOS spectrum. The emission 
from C2H2 at 13.7 /Ltm and HCN at 14 /itm is clearly no 
longer as strong as it was in quiescence. The case of 
CO2 is more hard to judge. The Q-branch of CO2 at 
14.97 nm. is close to a strong water line at ~14.9 /zm. 



and overlaps with one OH cross-ladder doublet (see Fig- 
and Table [2]) . It is therefore hard to say if a small 
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contribution from CO2 is still present in outburst, but 
the observations can be explained by H2O-I-OH emission 
alone. This is consistent with the non-det ection of CO2 



and C2H2 an d HCN as well) reported in Pontoppidan 



et al. (2010a 



for EX Lupi in outburst, using the C08 

spectrum. 

We then use our LTE model to add the organic 
spectral features in the outburst spectrum and derive 
upper limits. We start with the model parameters 
assumed in quiescence and decrease A^moi until the line 
fluxes get below the noise limit in outburst. The result 
is that all organic molecules have to be reduced in 
column density by at least a factor ~ 5 to be undetected 
in outburst. If the emitting area and/or the tempera- 
ture increase in outburst then this factor is ^ 10 or more. 



5.4. Summary of Changes in Emission from Model 
Results 

We report the constraints derived using single-slab 
LTE models in Table pi In Figures raandnil we show all 
together the LTE models derived tor water and OH and 
the models assumed for organics. By comparison of the 
model results in quiescence and outburst we can summa- 
rize the changes in molecular emission as follows. Water 
becomes more optically thin and the column density de- 
creases (iVn^o from « 3 X 10^^ to « 3 x 10^^ cm"^)^ 
but the emission comes from a larger area. Given the 
degeneracy between parameters, another possible solu- 
tion is that the temperature increases while the column 
density remains the same over the same emitting area as 
in quiescence. However, the larger emitting area of OH 
in outburst supports the former solution, assuming that 
the two molecules emit from the same por tion of the disk 
(as found e.g. by 
ing water and Ot 



Doppmann et al. 2011). So, consider- 
together we would conclude that the 
emitting area of warm gas increases from ~ 0.6 to ~ 1.3 
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TABLE 5 

Constraints on H2O and OH Emission from Single-slab LTE Models. 



Line ScimplG 




Parameter 


Quisscsnc© ^S05) 


Outburst (AOS) 




Tex (K) 


580+1?, 


595+"° 


H20 


A^mol (10^^ cm~2) 


32.0ti3;f 


2.8l?;? 






(AU) 


n c:o+0.07 
^■3»-0.08 


(1.30) 




Tex (K) 


(590) 


59011° 


low-Bu OH 




< 0.7 


2.8tS:^ 






(AU) 


(0.58) 


I.30IHI 



Note. — The confidence limits on the parameters are derived using constant boundaries. In particular we report the 
Ax^ ~ 1.0 boundaries, which in one-dimensional planes give the 68.3% confidence on each parameter estimate. Values in 
parentheses are assumed from water and OH emission that are well characterized (see Section Isl. For OH we consider here 
only lines with Eu < 6000 K (see Table [2]|. 

AU in radiufFM Remarkably, while the column density of 
water - and of C2H2, HCN and CO2 - decreases in out- 
burst, the OH column density increases. OH emission 
can be explained by two temperature components, one 
characterized by the excitation of very high rotational 
levels > 10, 000 K) and one by transitions with £"„ 
lower than ~ 6000 K. The second component shows a 
temperature that is consistent with the result found for 
H2O in outburst (Tex ~ 600 K) and increases in column 
density in outburst for at least a factor 3 (A^oh from 
< 0.7 X 10^^ to « 3 X 10^^ cm-2). 



6.1. H2O and OH Emission between Quiescence and 
Outburst 

During the 2008 outburst in EX Lupi the mass ac- 
cretion rate onto the star was found to increase from 
- 10-i°-10-9 to ~ 10"^ M0 yr"^ and the stel- 
lar -t-accretioiyuimnosi^ 

^ 4.5 ( |Sipos et al ||2009i lAspin et al.1|2010| [Juhasz et ah 
2011 1. Whatever be the heating source of the gas in the 
disk, water formation via gas-phase reactions proceeds 
vigorously when the gas temperature is above 3 00 K (e.g 



Glassgold et al. |2009[ [Bethell fc Bergin||2009| ) , and can 
be exp ected to increase durin g an outburst. According 



6. IMPLICATIONS 

Our study demonstrates that accretion outbursts in 
EX Lupi considerably affect the molecular chemistry 
in a region that is probably the surface of the disk 
at radii relevant for planet formation. If confirmed, a 
change in emitting area of OH and water would show 
that a larger extent of the disk is heated significantly 
during outburst. But a greater illumination and heating 
are not enough to explain the change in emission we 
observe, and we should consider the consequences of 
the spectral energy distribution of the emission. UV 
photodissociation of water is indeed the best candidate 
we find in the literature to explain the sudden formation 
of a large amount of OH in outburst. Moreover, if 
water becomes more optically thin in outburst and its 
column density decreases, this may allow UV radiation 
to penetrate deeper in the disk. Organics that might 
have been shielded by water in quiescence would not be 
so any more, and could be photodissociated. A thinner 
layer of water would also probably favor more emission 
from H2 within disk regions that were previously not 
heated during quiescence. In the following we briefly 
discuss the scenario we propose to explain the change 
in molecular emission observed in EX Lupi between 
quiescence and outburst, and we compare our results 
with theoretical predictions from the literature. 



to the Glassgold et al. (2009) model, the range of H2O 
column densities derived tor EX Lupi in the two phases 
(« 10^''-10^^ cm^^) can be explained with heating by 
X-rays irradiation alone, provided that H2 formation on 
grains is efhcient enough. However, it should be checked 
under which conditions formation of water would over- 
come its destruction under an increased X-ray radiation, 
and compare with the case of EX Lupi. Unfortunately, 
X-ray observations were not taken close enough in time to 
the Spitzer spectra we consider in this work, an d com- 
parison of the avai lable data from Gii del et al. (20101 
(August 2007) and [Grosso et al] ( |20T0| (August ^UUSJ 
does n ot provide conclusive constraints. [Glassgold et al. 
( 2009| also proposed that accretion (mechanical) heating 
would increase the H2O column density. This instead 
can be ruled out in EX Lupi: even if we had the same 
emitting area in quiescence and outburst the H2O col- 
umn density does not increase when the accretion rate is 
higher, more so if the emitting area is larger in outb urst 
as we assume in our reference case (see Section 5.1 ) 



11 Given the inclination of the disk of ~ 20° jSipos et al.|2009| 
these projected equivalent radii may be tracing actual physical disk 
radii. 



On the other hand, UV radiation has been proposed to 
play a key role in the formation/destruction and avail- 
abihty of H2O, OH, an d organic molecules in t he planet- 
forming zones of disks ( [Bethell &: Bergin|2009[ ). A strong 
UV radiation can photodissociate water vapor in favor 
of OH in the upper layers of the warm disk atmosphere, 
competing with the capability of H2O and OH to self- 
shield. The change in emission in EX Lupi is consis- 
tent with this scenario: we see N^^^q decreasing in out- 
burst, while A^oH increase up to the limit gi ven by the 
self-shielding mechanism (^ 2 x 10^^ cm~^, Bethell & 
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BerginpOOQ )■ Following in this direction, the increase of 
A'oH in outburst suggests that the connection between 
the strength of OH e mission and the streng th of the ac- 
cretion rate found by Carr & Najita (2011) may be at- 
tributed to the effect ot U V radiation rather than to me- 
chanical heating. Again, the available UV observations 
do not provide strong constraints but we have at least 
the evidence that UV radiation was higher at the time 
of the AOS spectrum than it was in quiescence, as in the 
U-ba nd (0.36 nm) the o bserved flux varied by a factor 
- 50 puhasz et al.pOll). 



The strongest evidence in favor of a key role of UV 
radiation probably comes from OH emission. The coex- 
istence of two temperature components in the MIR spec- 
trum of EX Lupi is an interesting possibility. In previous 
works using Spitzer spectra of classical T T auri systems 



a co ld component (T ^ 500 K) was found by Carr & Na- 
jita ( 2008j) fitting OH hues detected in the Lti module, 
while Carr & iNajita ( |2011 ) proposed a hot component 
(T ~ 4000 K) considering lines in th e SH module. Based 



on previous observational evidence, Najita et al.| ( |2010 1 
proposed that UV-irradiatcd disks may show a twofold 
OH emission: a prompt non-thermal component from 
UV photodissociation of H2O which would eventually re- 
lax to a thermal emission at the temperature of the gas. 
The broad population of OH le vels in the ra n ge Eu ~ 
400-28,000 K, first observed by |Tappe et ah] (j2008| in 
the HH211 outflow, was indeed explained with a radiative 
cascade down the ground-state rotational ladders follow 
ing UV photodissociation of H2O, which would pr eferen 



tially popu late OH levels with E„ > 40, 000 K ( |Tappe 
et al.||2008[ [Harich et al.||2000[ ). If our approximation ot 



two temperature components in outburst is confirmed, 
it would represent the first observational evidence of an 
ongoing prompt "hot" (in reality non-thermal) emission 
followed by a thermal emission at the temperature of the 
gas. In fact, in outburst we find one component that 
may hardly be explained with thermal excitation, while 
another component shows the same temperature found 
for water (see Section 5.2). An interesting peculiarity of 



EX Lupi is the detection of the cross-ladder transitions 
between ~ 15 and 19 ^m, that were never observed be- 
fore in any classical T Tauri system. The case of EX 
Lupi, as well as similar variable sources, might be very 
useful to clarify the nature of OH emission in the inner 
regions of circumstellar disks, but more theoretical work 
is needed. For instance, it would be interesting to know 
if the proposed excitation mechanism by radiative cas- 
cade would produce the curve that we see in the rotation 
diagram, and, if that is a transitory process, on which 
timescales/conditions it should be observed as two com- 
ponents rather than a continuum of temperatures. 

The simultaneous monitoring of UV, X-rays and MIR 
emission during future outbursts would help in confirm- 
ing the connection between the variation in accretion 
rate and changes in H2O and OH emission. Velocity- 
and transition-resolved spectral lines might be the key 
tool to clarify the location of the MIR molecular emis- 
sion as well as help in understanding the processes re- 
sponsible for formation/destruction of t he molecules we 



observe 



e.g. 



Pontoppidan et al. 2010b). More compre 



hensive chemical disk models accounting for both X-rays 
and UV variable irradiation would also be very helpful in 
better constraining water formation and survival in the 



inner zones of disks based on spectrally resolved synoptic 
observations of young variable star-|-disk systems. 

6.2. The Lack of C2H2, HCN and CO2 in Outburst 

The fact that organic emission dicreases in outburst is 
very interesting, and probably the most puzzling of our 
findings. Production of HCN has been proposed to be 
driven b y strong UV photodissociating N2 in disk atmo- 
spheres ( Pascucci et al. 2009) . In EX Lupi we observe 
instead a lower MCIN column density in outburst, and yet 
the ratio Nyich / Nq^yi^ ~ 2.7 in quiescen ce is consistent 
with t he median sun-like star reported in Pascucci et al. 
(l2009|. The lack of C2H2, HCN and CO2 in outburst 



seems to preferentially imply the large photodissociation 
of these molecules by UV, in contrast with the capability 
of self-shielding of H2O and OH. If the organic emis- 
sion probes the same disk radii as H2O and OH, then 
their non-detection would show that warm organics are 
located higher in the disk atmosphere than an optically 
thick layer that could shield them, wheth er it is made 
of du s t or H 2O a nd OH as suggested in [Pontoppidan| 
et al.| ( |2010a[ ) and |Bethell fc Bergin| ( |2009[ ) respectively. 
'I'he unclear behavior ot the emissi on from these organic 
molecules in other T Tauri disks (Carr & Najita 20111 
suggests that we still do not fully understand where ex- 
actly their emission originates and what are the factors 
most responsible for their abundance. It would be very 
interesting to see if organics are detected in EX Lupi 
again after the outburst. If that was found to be true 



than the vertical mixing proposed by Juhasz et ah] (2011 1 
to be powered- up in outburst may be the mechanism 
replenishing the disk atmosphere with molecules from 
closer to the midplane. Monitoring EX Lupi during the 
current new quiescent phase and future outbursts may be 
extremely important to check the timescales/processes 
relevant for organic chemistry in planet-forming regions 
of disks. 



6.3. Note on HI and H 



2 m 



Outburst 



We do not focus in this paper on the strong lines from H 
I an d H2 that appear in EX Lupi in outburst (see Section 
3.1 1. The increase in the accretion rate is certainly a 
possible explanation for an increase in H i emission (e.g. 
Pascucci et al.||2007 ), but photoevaporative winds might 



also c ontribute (see e.g. the discussion in 
2010 1. The increase in H2 emission can 



Najita et al 
3e explained 



by em ission from a warm disk atmosphere (Najita et al 



2010| iGorti et ar]|2011[ ), which in EX Lupi in outburst 
IS probably warmed up in regions that were too cold in 
quiescence (larger radii and/or deeper layers). It would 
be interesting to see whether a direct relation between 
these species and the formation/destruction of H2O in 
EX Lupi exists. If warm H2 is indeed available on a 
larger portion of the disk at mosphere, then water v apor 
formation might be favored (Glassgold et al. 20091 and 
increase again the observed column density after the end 
of the outburst, when UV photodissociation diminishes. 

7. CONCLUDING REMARKS 

The molecular emission we detect in the quiescent 
S05 spectrum (from H2O, OH, HCN, C2H2 and CO2) is 
comparable to what has been found to be common in 
several other T Tauri systems ( Pascucci et al.|2009} Carr| 
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& Najita 2008 2011). This is a further confirmation 



that EX Lupi IS not distinguisha ble from typ i cal T Tauri 



systems, as already suggested by Sipos et al. (20091 from 
consideration of the spectral energy distribution m qui- 
escence. Therefore, the remarkable changes in molecular 
emission observed in the outbursting EX Lupi might 
be extremely important to constrain chemical-physical 
processes common to all T Tauri systems. In this paper 
we have addressed the simple scenario of the increase of 
one system parameter in EX Lupi; the disk-atmosphere 
illumination during outburst. To first order this is con- 
sistent with an outbursting event that accreted material 
only from inside an inner hole in the dusty disk (~ 0.3 



AU) and not over global disk scales 


Juhasz et al. 


2011 


Goto et al. 


2011 Kospal et al. 2011 


. Considerat 


ion 



more sophisticated scenarios needs to be included as 
soon as high-resolution data will become available. The 
geometrical structure of the EX Lupi system is probably 
affected by strong outbursts and may be relevant for 
the correct interpretation of the changes in emission we 
observe, including the increase of line fluxes in the LH 
module. As the stellar -|-accretion luminosi ty L^. + Lgcc 
increased in outburst by a factor ^ 4.5 (Aspin et al 



20101, the dust-sublimation radius is probably shitted 
outward by a factor ~ 2. This can have effects on the 
inner rim location and the shadowed portion of the disk 
as well, and/or the location of the snow lines of the 
different chemical species. The location of the disk we 
likely probe, the inner few AUs, is indeed where both 
the shadowing and the snow line are believed to be. In 
addition to that, UV photodesorption from icy grains 
is expected to power the production of OH as well as 



to largely replenish water vapor in the disk atmosphere 
Despite the increased UV and the 



(Oberg et al. 


2009 


production o 


OH, 



of water vapor in EX Lupi in outburst. If self-shielding 
of water breaks under particularly h arsh conditions (see 
the cas e of Herbig AeBe stars, e.g. I Pontoppidan et a~ 
( 2010 a [) and Fedel e e t al.| (2011), and transitional disks, 
e.g. JNajita et al. (|2010 )), then a receding snow line 
may still provide the ice reservoir needed to sustain 
a wet disk atmosphere in EX Lupi. Velocity-resolved 
observations are critical to confirm the location of the 
different molecular emissions and explore these senarios. 
New infrared facilities from the ground (e.g. the soon 
upgraded VISIR on the Very Large Telescope, and in 
the future METIS on the E-ELT) and from space (with 
MIRI on the James Webb Space Telescope and MIRES 
on SPICA) will be the key tool in future investigations 
of the warm molecular emission from the inner regions 
of disks. 
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APPENDIX 
A SINGLE-SLAB LTE MODEL 

A single-slab model is used to extract from the spectra the excitation conditions of the emitting gas. We assume the 
gas to be in local thermodynamical equilibrium (LTE), but account for optical depth effects. We follow a formalism 
similar to van der Tak et al. (2007), and we ignore the absorption of background radiation by the emitting gas (i.e. 
the absorption term of the radiative transfer equation). Neglecting mutual overlap between transitions and assuming 
the line profile function to be Gaussian, the integrated intensity (e.g. erg cm~^ sr~^) of a transition connecting 
the molecular level u (upper) and I (lower) reads 



Td — 



Aw,. 



Bu{Tcyi) 



1 



exp 



dy. 



(Al) 



Here, r„; is the opacity at the line center, B^(Tcx) the Planck function for the excitation temperature Tqx, Vui the line 
frequency and Av the line width, which is fixed to 1 km s^^. The integral in the above expression has been obtained 
using a change in variable y = 2-^/ln( 2)f /At; and approaches \/7r ( 1 — e""^"') for TuI < 1. Thus, the usual expression 
used in escape probability codes (e.g. van der Tak et aL]|2007 ) is recovered for low opacities. The opacity at the line 
center is given by 

^-^"'V (.,'-^~..), (A2) 



Tul 



An^/Tr Avvf, 



91 



with the Einstein- A coefficient of the transition yl„; and the total molecular column density N^^oi. The statistical 
weight and the normalized level population (X^i^i = 1) of the upper and lower level are given by gi and x^, 
xi, respectively. The level population is obtained from the Boltzmann distribution Xi = gtexp (—Ei/kTcx) /Q(Tcx), 
with the partition sum Q (Tex) and the energy of the level Ei. The relevant molecular parameters are taken from the 
HITRAN 2008 database ( [Rothman eraI1|2009[ ) . 

To model the Spitzer spectra, i.J the intensity of all transitions of a chosen molecule (e.g. H2O) in the wavelength 
range of the Spitzer IRS is calculated using equations Al and A2 for an assumed Tgx and A^mob the integrated 
flux {Fui, in units of erg s~^ cm~^) is calculated by multiplying the intensity /„; with the solid angle AQ = A/(P, 



17 



where A is the emitting area in the sky, for which we assume a projected circle with radius Vp = y /A/n, and d the 



distan ce to the source (assumed to be consistent with the distance of the Lupus complex, 150 pc, [Lombardi et al 
2008| ; Hi.) a synthetic spectrum (in erg cm~^ fim"^) with spectral resolution R = A A/ A 600 is generated based 
on ail transitions in the wavelength range considered, with flux given by 



1 



(A3) 



where a\ — A/2i?-\/2 ln(2), and A^j and F^; are respectively the wavelength and the flux of a single transition. 
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